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ABSTRACT 

We investigate the energy release due to the large-scale structure formation and the subsequent transfer 
of energy from larger to smaller scales. We calculate the power spectra for the large-scale velocity field 
and show that the coupling of modes results in a transfer of power predominately from larger to smaller 
scales. We use the concept of cumulative energy for calculating which energy amount is deposited into 
the small scales during the cosmological structure evolution. To estimate the contribution due to the 
gravitational interaction only we perform our investigations by means of dark matter simulations. The 
global mean of the energy transfer increases with redshift ^ (z -I- 1)^; this can be traced back to the 
similar evolution of the merging rates of dark matter halos. 

The global mean energy transfer can be decomposed into its local contributions, which allows to 
determine the energy injection per mass into a local volume. The obtained energy injection rates are at 
least comparable with other energy sources driving the interstellar turbulence as, e.g. by the supernova 
kinetic feedback. 

On that basis we make the crude assumption that processes causing this energy transfer from large to 
small scales, e.g. the merging of halos, may contribute substantially to drive the ISM turbulence which 
may eventually result in star formation on much smaller scales. We propose that the ratio of the local 
energy injection rate to the energy already stored within small-scale motions is a rough measure for 
the probability of the local star formation efficiency applicable within cosmological large-scale n-body 
simulations. 

Subject headings: large-scale structure of universe, intergalactic medium, galaxies: interactions, stars: 
formation 
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1. INTRODUCTION 

During the last decade great success could be achieved 
in our understanding of the detailed mechanisms for the 
formation and evolution of cosmic structure. The role of 
the underlying cosmological models and parameters have 
been investigated by numerous numerical n-body simula- 
tions. On large scales the calculated distribution of mat- 
ter is in excellent agreement with the observed one, e.g. 
with the galaxy distribution and the distribution of the 
intergalactic medium (IGM). The fast enhancement of the 
available computational power permitted to cover an in- 
creasing dynamical range and/or to consider additional 
processes. In particular hydro dynamical models have been 
developed very successfully. 

The evolution of cosmic matter can roughly be subdi- 
vided in terms of scales: On spatial scales larger than the 
Jeans length gravitation dominates and on smaller scales 
hydrodynamical processes do. On the other hand for the 
description of the IGM, e.g., this distinction is insufficient. 
The IGM is strongly influenced by both, the large scale 
structure evolution and the feedback of the luminous mat- 
ter, in particular by the star formation processes: Super- 
nova explosions sweep out the galactic gas enriched by 
heavy elements into the IGM changing its chemical com- 
position and thermal state. Radiation ionizates the IGM 
in the environment of the galaxies. 

In order to obtain an appropriate description of the 
physical state of the large-scale distributed gas also the 
amount and the distribution of stars and their back-reaction 
has to be estimated. However, to link the process of star 
formation to the large-scale structure evolution is by far 
out of scope. 

Incorporating the stellar feed-back in simulations in- 
evitably needs to connect the star formation rate to avail- 
able gas parameters. Those parameters can be, e.g. the 
local density and the local gas temperature. Note, in terms 
of the considered simulations 'local' stands for the average 
over smallest resolved scales, usually of the order of 1 - 100 
kpc. Schmidt (1959) found that in the interstellar gas the 
star formation rate p* (SFR) is related to the density p by 
p* oc p", where n is adopted to be about 1.5 (cf. Kenni- 
cutt 1998, and references therein). Although those scales 
are far below the resolution in large-scale simulations, this 
empirical relation is often applied. Gas which fulfills cer- 
tain density and temperature criteria is assumed to form 
stars according to the above given Schmidt-law (e.g. Yepes 
et al. 1997; Springel 2000; Steinmetz 2001; Nagamine et al. 
2001; Ascasibar et al. 2001). An alternative approach 
has been introduced by Kauffmann et al. (1999) linking 
semi-analytic galaxy models to bound dark matter halos. 
Applying such prescriptions for star formation, e.g., per- 
mits to calculate the star formation history and the stellar 
metallicity distribution in the universe. However, the va- 
riety of used criteria (cf. Kay et al. 2001) indicates that 
the conjunction is still uncertain. 

Knowing the processes on the scale of Molecular Clouds 
(MCs) which most probably arc controlling the star for- 
mation rate would provide an indication on the possible 
linking quantities. Therefore, let us shortly summarize 
some recent results of detailed investigations related to 
star formation. 

Star formation is hosted by interstellar clouds of molecu- 



lar hydrogen. Stars probably arise from shock-compressed 
dense cores within the clouds (Blitz & Williams 1999). 
The cores are produced by supersonic motions of the gas 
due to the presence of turbulence (Burkert 2001). Accord- 
ing to Klcsscn (2000) the SFR especially depends on the 
scales on which turbulence is driven. 

There are indications that the formation of clouds is 
linked to larger scales: Blitz & Williams (1999) argued 
that MCs are formed through the condensation of H I re- 
gions in conjunction with some other mechanism as has 
been proposed by Ballesteros-Paredes et al. (1999), e.g. 
due to colliding interstellar gas streams (cf. Burkert 2001). 

The hfe-time of the MCs, 20 - 100 Myr, indicates that 
the internal turbulence is more or less continuously driven 
by external forces (Mac Low et al. 1998; Padoan et al. 
2001). Nevertheless, star formation occurs probably only 
once within one cloud (Elmegreen 2000). Consequently, 
the rate of star formation may depend on the presence of 
an effective intercloud turbulence. 

The state of the ISM, i.e. whether it is turbulent or not, 
may be derived from the power spectrum of the velocity 
field or from that of the matter distribution. For instance, 
Goldman (2000) investigated the neutral hydrogen distri- 
bution of the Small Magellanic Cloud (SMC) and derived 
the spatial power spectrum. He found that the spectrum 
of the SMC is to some extent steeper than Kolmogorov's 
1/3-law and explained that by the compressibility of the 
gas. Furthermore, he suggested that the origin of the tur- 
bulence could be caused by a close encounter between the 
Small and the Large Magellanic Cloud (LMC) 2 x 10^ yr 
ago and he showed that supernova feedback could only 
contribute little to the kinetic energy of the SMC. This 
illustrates one possible mechanism to force the ISM into a 
turbulent state: kinetic energy is transferred from larger to 
smaller scales, i.e., kinetic energy from the relative motion 
of the SMC and LMC is transferred into internal motions. 
It is intriguing in respect to the relation between turbu- 
lence and star formation that both the SMC and the LMC 
show evidence for large increases of their star formation 
rate at the time of the encounter (Larson 2001). 

From these investigations one may conclude that the 
SFR is controlled by the rate of formation of dense sub- 
structures - or clouds - in the interstellar medium, dense 
enough to undergo gravitational collapse and fragmenta- 
tion. Whether these substructures are generated by ther- 
mal instability, gravitational instability, turbulent com- 
pression or some other process is still under debate (cf. 
introduction in Scalo & Biswas 2001). We suppose that 
interstellar turbulence is the leading mechanism. Then 
the question arises, how is the interstellar medium forced 
into a turbulent state? What drives the turbulence over 
a sufficient long time period needed to form eventually 
stars? Again, different energy sources such as supernova 
feedback, galactic differential rotation (Sellwood & Bal- 
bus 1999), infaU of high- velocity clouds (Blitz et al. 1999) 
and other mechanism are conceivable. Merging events be- 
tween galaxies or minor-merging events between galaxies 
and smaller clouds are supposed to cause efficient star for- 
mation (Kolatt et al. 1999; Somerville et al. 2001; Kauff- 
mann et al. 2001). This may indicate that merging, which 
occurs due to the structure formation, could provide suf- 
ficient energy to stir the interstellar medium. 
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Therefore, the question arises how much energy is even- 
tually available by injection from the extra-galactic scales. 
When large-scale structures are formed gravitational en- 
ergy is released and stored in large-scale motions. The 
structure evolution transfers the kinetic energy from large- 
scales to galactic scales. May this transferred energy sig- 
nificantly contribute for balancing the dissipation of the 
turbulent field on the sub-galactic scales? 

In this paper we want to address the question which en- 
ergy transfer down to galactic scales can be expected from 
the large-scale structure formation, what is the time evolu- 
tion of the transfer and how is it spatially distributed. On 
large scales (> 0.5 Mpc) the baryonic matter is still tightly 
coupled to the dark matter, thus the energy transfer is gov- 
erned by the gravitational interaction only. We use cos- 
mological n-body simulations to determine the transfer. 
We start with the investigation of the power spectrum of 
the large-scale velocity field, closely related to the energy 
spectrum (Sec. 2). Then by introducing the concept of 
cumulative energy we determine the energy injection into 
galactic-scale motions (Sec. 3). We calculate the volume 
averaged energy transfer and discuss its relation to the 
merging rate evolution of halos (Sec. 5). Furthermore, the 
spatial distribution of the energy injection at given scale 
is determined. This allows to infer a local energy injection 
rate which can be attributed to a heating rate with respect 
to the baryonic matter (Sec. 6). 

Making the crude simplification that the obtained en- 
ergy injection rates at galactic scales are linked to the nec- 
essary energy input for driving the interstellar turbulence 
and moreover assiiming that this turbulence is controlling 
the star formation leads us to an estimator for the local - 
in terms of cosmological simulations - star formation rate 
(Sec. 7). Finally we summarize our results and discuss 
possible implications (Sec. 8). 

2. POWER SPECTRUM OF THE LARGE-SCALE VELOCITY 

FIELD 

We calculate the power spectrum of the large-scale ve- 
locity field using n-body simulations of the cosmic struc- 
ture formation. The spectrum should possess similar fea- 
tures as the power spectrum of the density field: At scales 
above « 60 Mpc (at present) the modes are decoupled and 
still evolve according to linear growth whereas at smaller 
scales the evolution of the modes is highly non-linear and 
coupled. 

We simulate the formation of the large-scale structure 
using a particle-particle/particle- mesh (P"^M) code consid- 
ering only the interaction of the dark matter. We have 
performed diS^erent simulations with box sizes L\)OK from. 
8 to 250 /i"^ Mpc and with 256^ and 128^ particles. We 
apply a cosmological model with matter density Qq = 0.3, 
cosmological constant ^Ia = 0.7, baryon density fis = 0.04 
and a Hubble constant Hq = hH^°^ = 70kms-^ Mpc"^ 
In this paper we are mainly interested in establishing a 
method to calculate the local energy release due the for- 
mation of the large-scale structure. We therefore apply 
only the customary ACDM model given above. 

Using n-body simulations the velocity field is only de- 
fined at particle positions. Therefore, calculating the power 
spectrum of the velocity field by its Fourier transform 
would require to interpolate the velocity field to a grid. 
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Fig. 1. — Power spectra of the velocity field at different redshifts 
z calculated by the algorithm introduced by Lomb (1976). The 
solid line gives the spectrum according to the linear growth at z = 
0.05 The vertical dotted line gives the Nyquist wavenumber which 
corresponds to the minimal distance of two particles if the particles 
are evenly distributed. 



We avoid this procedure by applying an algorithm which 
was introduced first by Lomb (1976) to calculate the 'pe- 
riodogram' of a randomly sampled time-dependent signal. 
We apply the algorithm in its form given by Press et al. 
(1986) and extend it to three-dimensional data sets. 

The most likely power P„ (fc) of a mode with wavevector 
k for a velocity field sampled at N points Xj {j e N) is 
given by 

2 

Ej(^j cosej(fe) 



1 



+ 



e,(fe) 



E,cos2e,(fe) 



J2Avj -v)smej{k) 



E,sin^e,(fe) 



(1) 



^ J. Y^i sin 2fca;i 

-arctan^z^^^ ^7— - 

2 2^ . cos 2kxi 



where Xg denotes the mean volume per sampling point. 
Averaging the power of several arbitrary wavevectors with 
the same absolute value |fc| = k eliminates the angular 
dependency 

P,(fc) = (P,(fe))|fe|^,. 

If the velocity field is evenly sampled the algorithm is iden- 
tical to the calculations of the power spectrum by the 
Fourier transform. If the field is unevenly sampled the 
algorithm is equivalent to a harmonic least-squares anal- 
ysis (Scargle 1982). Therefore, the largest wavenumber, 
power of which can unambiguously determined, largely ex- 
ceeds the Nyquist wavenumber of an evenly spaced parti- 
cle distribution. Figure 1 shows velocity spectra of one 
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simulation taken at different redshifts z. A Fourier trans- 
form based power spectrum always rapidly decreases when 
approaching the Nyquist wavenumber, whereas the algo- 
rithm of Lomb does not. The P^M simulations resolve ve- 
locity modes whose wavelength is smaller than the Nyquist 
one and the algorithm of Lomb allows to calculate its 
power. 

At redshift zero the spectra show at fc « 0.2 h Mpc~^ 
the transition from linear to nonlinear growth alike the 
spectra of the density field. Alodes with wavelength above 
this transition scale evolve independently and according 
to the linear growth function, whereas modes with smaller 
wavelength are coupled to each other and increase much 
faster. This becomes more conspicuous by combining the 
spectra of simulations for different box sizes, see Fig. 2. 
At large scales again the linear growth is reproduced ('lin- 
ear region'), i.e. the modes evolve independly. Towards 
smaller scales the slope decreases and a nearly 'flat' re- 
gion occurs. After the flat region the spectra probably 
exhibit a common, power law-like behavior ('asymptotic 
region'). This deviation from the linear growth in the 
flat and asymptotic region characterizes the coupled mode 
evolution. The power throughout the flat and asymptotic 
regions depends on the largest modes within the simula- 
tion box: Figure 2 shows that the power of the largest 
wavelength Amax realized in a simulation evolves always 
according to the linear growth, indicated by the solid line. 
The largest wavelengths of the depicted power spectra are 
located either in the linear or in the flat region. The two 
spectra with Amax in the linear region exhibit roughly the 
same power in the flat region. The two spectra with Amax 
in the flat region exhibit in contrast a much smaller power 
throughout the flat and asymptotic region; the amplitude 
is governed by the linear growth of the power at the largest 
wavelength Amax- This shows that in the asymptotic and 
flat regime the power of a given mode is dominated by the 
larger modes. There is no evidence that smaller modes 
are able to increase the power of the larger ones, i.e. the 
transfer of energy is one-directional from larger to smaller 
scales. 

In the asymptotic region the shape of the spectra can 
be described by a power law with an index of about -4 
to -4.5. The overlap of the different spectra in this re- 
gion reflects an asymptotic behavior which is independent 
from the largest modes in the box. Qualitatively the same 
three regions can be distinguished in the power spectrum 
of the density field using high resolution simulations (cf. 
the spectrum given by Klypin et al. 2000). 

The power spectrum P„(fc) of the large-scale velocity 
field can be transposed into a velocity-scale relation by 
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(2) 



The asymptotic region with power-law behavior oc k'^ re- 
sults in the relation v oc with k = — (3 -I- x)/2 « 0.6. 

The asymptotic power-law behavior indicates on self- 
similarity at small scales, which should result in a scaling 
relation, connecting scale velocity vi, and energy transfer 
e. We provide a heuristic derivation of a scaling relation 
based on dimensional analysis. 

The energy rate per volume element and time has the 
following unit: [e] = erg cm~'^ s~^ = g cm~^ s~^. The 
driving force for all processes described by the simulations 
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Fk;. 2. - - Power spectra of the velocity field taken from simulations 
with different box sizes at same redshift z = 0.33. The labels denote 
box sizes in [/i^-'Mpc]. The simulation with label (L) was calculated 
with 256"* particles; the other simulations were calculated with 128"^ 
particles. The solid line shows the power according to the linear 
growth. The dotted line gives an estimated asymptote. 



is gravitation and, therefore, the gravitational constant G 
must enter the expression for the energy transfer e. Fur- 
ther quantities characterizing the behavior of the n-body 
system is the velocity vi, the associated length scale I, 
and the density p. The only possible combination of these 
quantities including the gravitational constant G to build 
a quantity with units [e] is 

' (3) 



e oc 



Gl^' 



The assumption of a scale independent transfer rate leads 

to the relation between the velocity dispersion and the as- 
sociated length scale vi oc l^^^, which matches the asymp- 
totic power law obtained by the simulations. The steeper 
decline of the scaling-relation, compared to Kolmogorov's 
relation vi oc may be attributed to the gravitational 
interaction leading to a contribution of compressional modes. 

The above results support the idea of a nearly constant 
energy transfer throughout the scales of the asymptotic re- 
gion. We wish, however, to obtain an exact determination 
for the energy transfer rate through a given scale. This 
will be provided in the next section. It is worth to men- 
tion some restriction of Lomb's algorithm if applying to the 
power spectrum of the velocity field in n-body simulations: 
Lomb's method rests on the assumption that the sampling 
points are randomly distributed, whereas the particles in 
the simulation are mainly located in the dense regions. 
Therefore, the true spectrum might slightly differ from 
the spectrum obtained by using Lomb's method. However, 
since the velocities are weighted by the local density the 
obtained spectra reflect rather the correct power spectrum 
for the kinetic energy. In this respect they where used and 
interpreted, indeed. A more detailed investigation of the 
velocity power spectra for different initial conditions and 
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cosmological parameters will be given in a forthcoming 
paper. 

In the following we use two results of the investigation of 
the velocity power spectrum: First, the large-scale velocity 
modes dominate the power of modes at smaller scales and, 
second, when forming the large-scale structure the energy 
transfer is mainly directed from larger to smaller scales. 

3. SCALE-BY-SCALE ENERGY BUDGET 

We want to calculate the gain of energy stored in mo- 
tions with wavelengths below a given cut-off \k = 2tt/K. 
For that purpose we utilize the concept of cumulative en- 
ergy which was introduced by Obukhov (1941), for a re- 
view see Frisch (1994). The cut-off wavenumber K divides 
the Fourier decomposition of the velocity field into a low- 
pass and a high-pass filtered part, consequently the field 
can be divided into a low-pass and a high-pass filtered field 



v{x) 



■■ — / dkv 
■■v{x)<+v{x)>. 



\k\>K 



dkv{k) 



(4) 

If V represents the velocity field of an incompressible medium 
the cumulative energy stored in the high-pass filtered ve- 
locity modes is given by 



C-Xin) 
■'K 



(5) 



where ( • ) denotes the volume average. Cumulative energy 
£^ can be gained or lost by energy injection by force J-^ 



K 



or by energy exchange 11^ between the high-pass and the 
low-pass filtered part: 



dtS 



>(in) _ -r> TT> 



(6) 

Dark matter does not exhibit any dissipation just by defini- 
tion but interacts via gravitation only. Thus applying the 
concept of cumulative energy to the density and velocity 
field of the dark matter the scale-by-scale energy budget 
Eq. (6) allows to calculate the energy exchange between 
different modes and, in particular, the energy transfer to 
the small-scale spectral range. 

The energy exchange 11^ gives an appropriate prescrip- 
tion how to determine the energy transfer: We have shown 
in the preceding section that the energy stored in motions 
in the 'asymptotic' region is governed by the transfer from 
larger to smaller modes, consequently we neglect the en- 
ergy input by force in the case of the high-pass filtered 
energy The high-pass filtered energy is only changed 
by energy exchange between the high-pass and the low- 
pass filtered part. To determine the energy exchange 11^ 
wc must calculate the time derivative of the high-pass fil- 
tered energy. 

Due to the properties of the cosmic matter the above de- 
scribed concept must be generalized to the case of a com- 
pressible medium. Thus the definition of the cumulative 
energy must be changed to allow also for the variation of 
the density field. Since the volume average ( • ) must now 
be applied to the product of three fields pv-v different pos- 
sibilities exist formally to construct the filtered energy. In 
the case of incompressible matter all possibilities of arrang- 
ing the filtering procedure result in the same cumulative 
energy = {v- v^) = {v^- v^) = {v^-v). Determining 
the cumulative energy of a compressible medium we have 



to generalize the definition Eq. (5), e.g. by {pv^ ■ v^) or 
{Pk '^k ' '^k)- show that the definition 

£^ = {iPv)^c-v>K) (7) 
is free of contradictions and we apply this to calculate the 
energy exchange. 

Wc evidence that the definition Eq. (7) provides a rea- 
sonable generalization on the basis of the dynamical equa- 
tions of the cosmic matter. We obtain well defined expres- 
sions for both the energy injection by force and the energy 
exchange. This is not possible for the other combinations 
of filtered velocity and density. The energy injection and 
exchange can be expressed by the filtered fields without 
time derivatives; to obtain these expressions we start with 
Eq. (6), insert the new definition for the cumulative energy, 
eliminate the time derivatives by the dynamical equations, 
and sort the resulting terms. 

The dynamics of the dark matter can be described by 
the Euler equation, i.e. the Navier-Stokes equation with- 
out viscosity, and the continuity equation. Introducing 
comoving coordinates distances are scaled by the Hubble 
expansion r = ax and peculiar velocities u are introduced 
hy V = Hr + au, where a denotes the scale factor which 
obeys the Friedman equation and which is normalized to 
unity today. Furthermore we introduce the 'comoving' 
density pc = p/a^ and eventually we get the dynamical 
equations 



dtu + (u ■ V)u + 2Hu = -^V$ 
dtPc + V • {pcu} = 0, 



(8) 
(9) 

where $ denotes the gravitational potential. Executing 
the steps given above we obtain the cumulative energy 
injection by force and the energy exchange in terms of the 
filtered fields 



- n< • 



-^((Pc«^ + r7^)- V$> 

( {pcU^ + {PcU)k) ■ {u ■ V)m ) + {u^- ^(V • PcU) ) 

In addition we introduced the time derivative operator 
Dt = dt + 411. The term AH takes into account the time 
dependence of the peculiar velocities due to the Hubble 
expansions and avoids the occurrence of a fictitious en- 
ergy exchange. Thus, we rewrite the scale- by- scale energy 
budget in the case of the cosmic matter 

(10) 

the individual terms are given above. 

We specify the demand for a 'well defined energy in- 
jection and exchange' by giving some requirements which 
have to be fulfilled. 

i. Additivity of the cumulative energy £ = £^ + £^ 

ii. Additivity of the cumulative injected energy by force 

J- = J-^ + 

iii. For the energy flux must be valid H^ + H^ = 

iv. Assuming p = const, the relations must go over into 
the incompressible case 

Using the definition of the cumulative energy Eq. (7) and 
the above derived expressions for the energy injection and 
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exchange it can straightforwardly be proven that the re- 
quirements i. to iv. are fulfiUed. Defining in contrast the 
cumulative energy by another combination of filtered fields 
would not fulfill these relations. 

Wc proceed with calculating the energy exchange using 
the high-pass filtered scale-by-scale energy budget Eq. (10). 
As discussed above the high-pass filtered energy injection 
by force can be neglected, thus we can calculate the 
energy exchange by 

n^ = -n> =A^^ = iA (11) 

This allows to determine at each moment the gain of en- 
ergy stored in small-scale motions when simulating the 
formation of the large-scale structure. We distinguish be- 
tween the energy exchange 11^ and the energy transfer 
rate e introduced above. The latter is usually introduced 
when treating an isotropic, quasi-stationary process of en- 
ergy transport and is the appropriate quantity if assuming 
a scale-independent energy transfer which allows to derive 
the velocity-scale relation. 

Contrary to turbulent baryonic matter, which possesses 
viscosity leading to dissipation at some scales, the behav- 
ior of the dark matter during the structure formation is 
completely different due to the absence of any dissipation. 
Once the energy is transferred to small scales it will not be 
dissipated but remains preserved as kinetic energy. Fig- 
ure 1 shows the resulting continuous increase of power in 
the velocity field. The energy exchange depends via the 
cumulative energy £^ on the cut-off wavelength \k, i.e. 
on the considered scale-length. 

As shown in the previous section the energy exchange 
is directed from large to small scales. Thus the energy 
exchange expressed by the cumulative energy delivers the 
energy gain of small-scale motions during the structure 
formation process. 

4. ENERGY EXCHANGE IN N-BODY SIMULATIONS 

To determine the cumulative energy in a n-body sim- 
ulation, which deals with a finite number of the particles, 
we rewrite the definition Eq. (7) in a discrete form. Each 
particle samples at its position a;„ the underlying fields, 
namely the velocity field M(a;„), the density of momentum 
pu{xn), etc.. We take into account that the fraction Vn of 
the integration volume, which is assigned to the particle 
n, depends on its distances to neighbored particles. We 
transpose the volume average ( • ) into a particle average 

We estimate the volume Vn by assuming an equidistant 
grid and dividing the cell volume VceW by the number of 
particles rii per cell i: Vn = Vce\i/ni. 

Using the additivity of the filtered fields, cf. Eq. (4), the 
high-pass filtered velocity field is given by = u — u^. 
The low-pass filtered field is obtained by calculating the 
average velocity within each cell u'j^{xi) = (tt)ccii, like- 
wise the density of momentum. The average field at the 
knots of the grid can only reasonably be determined if at 
least a few particles are within a given cell. Determin- 
ing the low-pass filtered field therefore connotes to intro- 
duce a density threshold. We use throughout this paper 
a low threshold, namely Ui > 3. In contrast to the veloc- 
ity field the density is only defined at the grid, therefore. 
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Fig. 3. — The evolution of the global energy exchange obtained 
from several simulations with different box parameters (box size and 
resolution, see caption Fig. 2) is shown. Due to the construction of 
the cumulative energy (cf. Eqs. (4) and (7)) the computed quantities 
depend on the cut-off wavelength \k- The latter depends in turn 
on the realized box size, since the mean particle distance is chosen 
as \fc. Throughout the asymptotic region the velocity power spec- 
trum behaves like Pu.{k) ~ A;~*'^^ which leads to the velocity-scale 
relation vi ~ l^^^ , (cf. Eq (2)). Since the high-pass filtered cumu- 
lative energy is roughly proportional to v^ the overall amplitude of 

the energy exchange scales approximately as Z^/^. This dependency 
can clearly be noticed within the shown figure. 



we can rewrite the low-pass filtered density of momentum 
(ptx)^(x,) = pQUi {u)ce\\{xi), where po denotes the back- 
ground density. Interpolating the filtered velocity field 
(u)ccii to the positions of the particles we can assign an 
'individual' dispersion to each particle 

al = {u{Xn) - (M)ceu(a;n)}^- (13) 

Inserting the definition of the individual volume, rewriting 

the filtered density of momentum, and using the dispersion 
per particle the definition of the discrete cumulative energy 
Eq. (12) leads to 

^'K = ~P.Y.yn, (14) 

where N denotes the number of cells. Thus, in n-body 
simulations the cumulative energy is easily obtained from 
the individual dispersions of the particles. We calculate 
the energy exchange, using Eq. (11) and the discrete form 
of the cumulative energy, when running the simulations 
introduced in Sec. 2. 

5. ENERGY EXCHANGE AND THE EVOLUTION OF THE 
MERGING RATE 

The energy exchange XV^ shows a continuous steep de- 
cline with time, see Fig. 3. The steep increase at the be- 
ginning of the 250 Mpc simulation is caused by the 
necessary density threshold to calculate the filtered fields. 
Therefore we do not discuss further this rather artificial ef- 
fect due to the finite mass resolution. Though the decrease 
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according to a power law witli the exponent 3 looks like an 
expansion effect, this is not the case. By construction, the 
energy exchange contains only the transfer between dif- 
ferent modes. Consequently Fig. 3 shows that the mean 
overall rate of energy exchange decreases with time. 

The energy exchange depends on the box size of the 
simulation, too, see Fig. 3. This is due to the fact that 
we choose the 'particle grid size' as cut-off wavelength. 
Increasing the box size entails that modes, belonging for- 
merly to the low-pass filtered modes, bcx;ome now part of 
the high-pass filtered ones. The differences in amplitudes 
can easily be estimated by using the asymptote at high 
wave numbers obtained in Sec. 2 and performing integra- 
tion with the corresponding cut-ofF length. The energy 
exchange from larger scales to this modes is additionally 
considered and increases the total amount of energy ex- 
change. 

To interpret the decreasing energy exchange we consider 
the common picture of hierarchical structure formation in 
more detail. At high redshifts small dark matter halos are 
formed the mass of which is continuously increasing by ac- 
cretion of 'free' matter and merging with other halos. Both 
processes lead to the transfer of energy formerly stored in 
large-scale motions into internal halo motions. For the 
time being we assume that the mean energy exchange per 
accretion and merging process is constant with time. Then 
the energy exchange rate reflects the number of accretion 
and merging events per time. Assuming furtliermore that 
the merger events are the dominant energy transfer pro- 
cesses, the evolution of the energy exchange ought to be 
equivalent to the merging history. The observed merging 
rates (Carlberg et al. 1994) and the rates determined in 
simulations (Gottlober et al. 2000) show indeed a behav- 
ior of ~ {z + 1)^. Thus the decrease of the mean energy 
exchange with time may be mainly attributed to the fact 
that merging of halos becomes more and more rare dur- 
ing the late structure formation. In the next section we 
will give further evidence that the assumption of approxi- 
mately constant mean energy exchange per accretion and 
merging event with time is reasonable. 

It is worth mentioning that in contrast to investigations 
elsewhere we are not forced to introduce a definition for the 
individual objects. The energy exchange traces the merg- 
ing of halos including a wide spectrum of object masses 
without the difficulty of identifying the individual objects 
and distinguishing between mergers and minor-mergers. 

6. SPECIFIC ENERGY EXCHANGE 

The expression for the global cumulative energy can be 
represented as the sum of local contributions: By rewriting 
its discrete definition Eq. (12) the sum over all particles 
can be decomposed into the sum over all cell contribu- 
tions each of which collects the contribution of all parti- 
cles within a given cell volume. That is possible because 
the cumulative energy was constructed by the high pass 
filtered quantities, i.e. by the local velocity dispersions. 
Then a local high-pass filtered energy density can be de- 
fined as 

where Xi denotes the position of the cell i and the sum is 
performed over all particles within the cell n G The con- 
tributions represent the energy content stored in high-pass 
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Fig. 4. — Volume averaged specific energy exchange taken from 
simulations with different box sizes (see caption Fig. 2). 

filtered motions with respect to the local volume related 
to the cut-off length. From the local energy density we 
can define the local contribution of energy exchange 

n>=\{D, + {u-V)}e>. (16) 

An example for the local energy transfer is the merging 
process of halos: When the separated halos move towards 
each other both of them carry kinetic energy attributed 
to large-scale motion. During the merging process at least 
part of this energy is transferred to the internal energy of 
the resulting halo. 

We introduce now the specific energy exchange F by 
dividing the local energy exchange ^^{xi) by the amount 
of mass within the volume of cell i. The specific energy 
exchange could be formally written as a heating rate with 
respect to the baryonic mass fraction. 



where nn denotes the number density of hydrogen atoms. 
Of course, this 'heating rate' cannot be directly interpreted 
as a true heating of the involved interstellar gas but serves 
solely as a convenient quantity for further estimates and 
comparisons. 

Figure 4 shows the evolution of mean the specific en- 
ergy exchange. This quantity depends also on the box 
size of the simulation for reasons which have been already 
discussed hitherto. The %h^^ Mpc simulation is missing 
power on large-scale modes, whereas the 250 Mpc does 
not resolve small halos. Therefore the remaining three 
simulations provide the best estimate of the specific en- 
ergy exchange. The evolution of the mean specific energy 
exchange F varies only slightly with time, compared to the 
behavior of the energy exchange itself. This supports the 
assumption that the typical energy exchange per merger is 
roughly c;onstant during the evolution of the universe and 
that the evolution of both the mean energy exchange and 
the merging rate are equivalent. 
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Fig. 5. Contours in density contrast p/po and specific energy exchange. Left panels are taken from the 80 /i ^ Mpc simulation at different 
redshifts. Right panels are tal^en from three different simulations at the same redshift z = 0.33. 
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Basically, the energy exchange only happens when ha- 
les merge or matter is accreted. Therefore, a correlation 
is expected between energy exchange and density. Fig- 
ure 5 shows that the specific energy exchange increases 
with density; the expectation value of the specific energy 
exchange is roughly proportional to the density. The en- 
ergy exchange is concentrated to a small fraction of the 
total volume; it is much more concentrated than the mat- 
ter, see Fig. 6. Since we consider the energy exchange per 
mass the description is likewise valid for the baryonic mat- 
ter as long it evolves coupled to the dark matter. This is 
also true for the motion of halos because the baryons are 
trapped in their gravitational wells. Thus the occurrence 
of a merging event and the resulting energy release can be 
well estimated by considering solely the dark matter. 

In contrast to the dark matter the baryonic one does not 
store all energy within small-scale motions. Instead, part 
of the energy injection at a given scale can be released as 
heat. The efficiency of that heating depends on the hy- 
drodynamical processes within the baryonic matter, i.e. 
mainly on the collision rate between particles. High den- 
sity acts in favor of both large stored kinetic energy and 
more effective dissipation into heat due to higher proba- 
bility of collisions. The partitioning of the injected energy 
into heat and kinetic motions is governed by the ratio of 
the infall velocity to the local sound velocity. Only if this 
ratio is sufficiently small, perturbations in the velocity field 
remain small and can be dissipated fast enough. 

The specific energy exchange provides the largest possi- 
ble amount of energy input into the baryonic matter aver- 
aged over galactic scale due to matter infall, merging, etc.. 
Figure 7 shows that the expectation value of the specific 
energy transfer rate F amounts to ss 10~^^ erg s"-'^ cm""^ 
when regions with density p/po > 300 are considered. The 
expectation value is proportional to the density, cf. Fig. 5. 
The density contrast within the galactic halo is at least 
two orders above 300, leading to a transfer rate about 



erg s ^ cm ^. 



10-24 

We want to estimate whether this energy injection could 
play a major role in the energy balance of the entire baryon 
content within a galaxy or not. The calculated energy in- 
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Fig. 7. — Distribution of the specific energy exchange into regions 
with density p/po > 300. 



jection rate is mainly caused by accretion and merging and 
to some extend by tidal action of neighbored halos. We 
wish to compare this released energy with the supernova 
(SN) feedback averaged over the same considered volume. 

From the above obtained specific energy rate the total 
energy injection rate can be obtained for a given volume 
containing the mass of a galaxy. Indeed, for a galaxy 
mass comparable with that of the Milky Way, i.e. of 
about 10"^" Mq, the total energy injection into the corre- 
sponding galaxy volume can be estimated by V^aiax F riH ~ 
V^gaiaxF (Mgalax/"^H/V'galax) = F Mgaiax/mR « 4 X 105°ergyr-^ 
The mechanical energy feedback per supernova is of the or- 
der of e X 10^^ erg. The efficiency parameter e is highly 
uncertain but surely below 0.5; most probably e is even 
much less (cf. Sellwood & Balbus 1999). Assuming a SN 
rate of 0.05 yr'^ (Timmes et al. 1997, found for the Milky 
Way a SN rate of 0.034±0.028yr~^) we get a mean energy 
injection rate of e x 10^^ ergyr^^. 

The above provided estimates lead to the conclusion 
that the energy input rate from large-scale evolution uni- 
formly distributed within a halo volume comparable with 
galactic size is certainly comparable at least with the en- 
ergy input rate by SN. 

7. ANOTHER PRESCRIPTION FOR STAR FORMATION 

For an appropriate description of the IGM in cosmolog- 
ical simulations the feedback of stars must be necessarily 
taken into account. Therefore, a prescription for the star 
formation rate per volume within simulations of the large- 
scale structure formation is needed. 

As briefly reviewed in the introduction it is usually as- 
sumed that baryonic matter which fulfills certain density 
and temperature criteria is able to form stars according to 
the Schmidt-law being aware that this is an empirical law, 
based on much different conditions. 

The 'true' conjunction between the available parameters 
in such simulations and star formation is still unknown 
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Fig. 8. — The global star formation rate according to Eq. (18), 
taken from simulations with box size Lbox = 80 h"^ Mpc (solid 
line). Depicted data points (cf. the compilation by Ascasibar ct al. 
2001) are from: PascarcUc ct al. (1998, p), Stcidcl ct al. (1999, m), 
Treyer et al. (1998, t), Lilly ct al. (1996, 1), Cowic ct al. (1999, c), 
Connolly et al. (1997, n), Sawicki et al. (1997, w), Gallcgo ct al. 
(1995, g), Tresse & Maddox (1998, d), Glazcbrook ct al. (1999, g), 
Yan ct al. (1999, y), Flores et al. (1999, f), Hughes ct al. (1998, 
h), Pascual ct al. (2002, u), Tresse ct al. (2002, c), Hopkins ct al. 
(2000, k), Sullivan ct al. (2000, v), Gronwall (1999, x), Moorwood 
ct al. (2000, o). Hammer ct al. (1997, r), and Haarsma ct al. (2000, 
a). We applied the correction given by Ascasibar ct al. (2001) to 
obtain star formation rates in the ACDM model. The solid gray 
line indicates the mean of the observed star formation rate. A quite 
different evolution is indicated by the data points given by Lanzetta 
et al. (2001) (hashed area). 



and allows to speculate and to set np very different mod- 
els. From detailed star formation models in turbulent MCs 
indications are provided that the star formation rate may 
be controlled by mechanisms which force interstellar tur- 
bulence. One possibility to stir the ISM may be the merg- 
ing of galaxies or the occurrence of minor- mergers. Such 
merging events are supposed to generate efficient star for- 
mation (e.g. Kauffmann et al. 2001). 

We have investigated the energy-injection into the inter- 
nal motions of the halos at sizes of « 30 kpc - 2 Mpc due 
to the large-scale structure formation and we have found 
that the energy per mass input is comparable with energy 
injection rates attributed to supernova. However, even if 
the energy transfer from the large-scale structure forma- 
tion provides sufficient energy to stir the ISM, at present 
it cannot be shown neither a mechanism to propagate the 
injected energy down to scales of observed ISM turbulence 
nor the detailed structure of the necessary turbulent field 
(for addressing this problem see however Jog & Ostriker 
1988; Sellwood & Balbus 1999). Although the effect of the 
energy transfer onto the ISM is unknown, it represents one 
of the possible sources for driving the interstellar turbu- 
lence. Therefore, we assume here that the energy input is 
the linking parameter and make the following proposition: 



The energy transfer from the large-scale strticture for- 
mation provides sufficient energy to stir the ISM and the 
transfer is the main energy source for interstellar turbu- 
lence. Then, as a consequence, the amount of star for- 
mation should be determined by this energy transfer as 
well. 

Following this quite speculative conjunction we would 
like to construct a prescription for the star formation rate 
within a small volume still described within large-scale 
structure formation simulations, i.e. related to scales of a 
few kpc. We assume direct proportionality to the local en- 
ergy injection rate ir^{xi) oc dt cr^- For a given halo 
part of the kinetic energy of the gas is transformed into 
heat. High temperature is expected to suppress star for- 
mation. We make the crude assumption that the thermal 
energy is a constant fraction proportional to the kinetic 
energy of the halo gas. This leads us to assume inverse 
proportionality of the SFR to the energy already stored 
within the small-scale motions e^{xi) cx X^nGi '^n- 

In the result we get for the fraction of mass which is 
transferred into stars during a certain time interval At 
within a given cell Xi 

7 — ^ = ^0 ^7 — r^i 

oc^^nei^At=^^^4^^t (18) 

where eg is a normalization constant which has to be de- 
termined, e.g. by the observed SFR at any fixed redshift, 
and m{xi) denotes the mass within the considered cell. As 
can be seen from the last term in Eq. (18) the ratio /e^ 
is independent from the local density. The comoving cos- 
mic SFR p» is obtained by summing up the contributions 
of all cells and dividing by the comoving volume of the 
simulation box L^^^^ 

Note, assuming for the SFR proportionality to the lo- 
cal specific energy injection the expression (18) represents 
the simplest combination of quantities obeying the dimen- 
sional requirements. 

Since the energy injection rate is closely correlated with 
the merging and the matter infall rates as was shown in an 
earlier section the star formation history calculated in this 
manner reffects mainly the inffuence of merging processes. 

For transparency we give the time evolution of the dif- 
ferent constituents of the expression (18) averaged over 
the comoving box volume, see Fig. 9. The amount of mass 
concentrated in regions with density contrast p/ pa > 3 in- 
creases continuously with time. The hierarchical cluster- 
ing process leads to increasing kinetic energy of the halo 
matter which happens according to the mass inflow and 
concentration. This explains the quite similar slope for 
the density and kinetic energy evolution. The energy in- 
jection into small scales is a slightly decreasing function 
with time as shown already in the preceding section, cf. 
Fig. 4. 

Implementing the prescription Eq. (18) within our simu- 
lations we obtain according to Eq. (19) the evolution of the 
cosmic star formation rate p*, as shown in Fig. 8. For com- 
parison we included a compilation of observational data of 
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star formation rates recently obtained for different red- 
shifts. Wc have fixed the normalization constant eg using 
the data in the vicinity of z — 0.5. Even if there is a 
big scatter in the estimated star formation rates it is com- 
monly accepted that there is a steep increase from redshift 
zero to redshift one (cf. gray line). Estimates at higher 
redshift indicate a subsequent decline even if these mea- 
surements are still under debate. Although dust extinction 
is to some extend allowed for it may still cause the decline 
of the data points. A quite different star formation history 
is recently given by Lanzetta et al. (2001). They found a 
continuous increase of the SFR with redshift. 

The evolution of the calculated 'star formation rate' 
matches the cosmic star formation history to some ex- 
tent. However, the limitations of the proposed model arc 
obvious: In addition to the considered impact of energy 
transfer, which mainly reflects the impact of merging and 
accretion, the star formation is affected by hydrodynami- 
cal effects, radiation, etc. related to the properties of the 
baryonic matter. In particular, the back-reaction of the al- 
ready formed stars onto the gas by winds, supernova etc. 
may locally turn off any subsequent star formation. These 
processes result almost exclusively in an reduction of star 
formation which may cause in particular a steeper slope 
at small redshifts than obtained here. 

The advantage of the proposed quantity is obvious as 
well: It can be computed on-the-flight without determin- 
ing the single matter halos and without tracing back their 
evolution. It can be used complementary within the frame- 
work of semi-analytical descriptions of star formation ap- 
plied to the results of n-body simulations. Of course, this 
cannot replace the knowledge of the detailed star forma- 
tion process including the full cooling history and detailed 
velocity fields. However, for the consideration of large- 
scale structure formation and the large scale distribution 
of quantities related to star formation this could serve as 
a suitable approximation. 



8. SUMMARY 

During the formation of the large-scale structures grav- 
itational energy is transferred into the internal motion of 
the dark matter halos. In particular by merging of the ha- 
los part of the kinetic energy is subsequently transferred 
from the large-scale movement into small-scale internal 
motions. The baryonic matter possesses the same amount 
of kinetic energy per mass as long as it moves togcthcir with 
the dark matter, i.e. as long as gravitational interaction 
dominates over pressure forces. Thus the large-scale struc- 
ture formation is transferring energy into small scales. We 
have addressed the question how this transfer can be de- 
scribed, how it evolves during the cosmological structure 
formation and whether this transferred energy rate is of 
comparable order of magnitude to compete with internal 
energy sources in the galaxies. 

We have performed cosmological n-body simulations to 
determine which energy amount the large-scale structure 
formation provides to a given small scale. 

A crucial point in that picture is how the energy transfer 
is directed. Wc have considered this question investigating 
the evolution of the power spectrum of the large-scale ve- 
locity field. To this end we have used a method introduced 
by Lomb (1976). The power spectrum evolves according 
to the linear growth at scales larger « 60 Mpc. Be- 
low these scales it possesses a more shallow region and at 
scales below « 0.5 Mpc it roughly behaves like a power 
law ~ In general, the features of the velocity spec- 

trum are similar to those of the power spectrum of the 
density field. The behavior of the velocity spectra shows 
that large-scale modes dominate the power of small-scale 
modes and that the energy transfer throughout the modes 
is almost exclusively one-directional, namely from larger 
to smaller scales. 

To calculate the energy transfer through a given scale we 
have used the concept of cumulative energy. The Fourier 
transforms of the velocity field and of the density of mo- 
mentum are subdivided into a high-pass and a low-pass 
filtered part. As a result the energy density can also be de- 
composed into a high-pass and a low-pass filtered part and 
the scale-by-scale energy budget contains a term describ- 
ing the energy exchange between the two spectral ranges. 
The exchange rate may be calculated by help of the time 
derivative of the high-pass filtered energy and equals to 
the energy injection rate into scales smaller than the given 
cut-off length. 

Using the involved quantities we can define straightfor- 
wardly the mean global energy exchange rate into given 
scales at any time during the performed n-body simula- 
tions. We obtain the mean energy exchange rate which in- 
creases with redshift according to a power law ~ (z -I- 1)^. 
We argue that this can be attributed to the mass increase 
of dark matter halos. The evolution includes all kinds of 
infall, i.e. continuous accretion of matter onto halos, minor 
mergers, or even merging of halos. If all processes evolve 
similar or if merging is dominant the evolution of the mean 
energy exchange exhibits the same time-dependence as the 
merging rate of dark matter halos. 

Due to the additivity of the cumulative energy with re- 
spect to its local contributions if applied to grid based 
n-body systems a local high-pass filtered mean energy ex- 
change rate can be defined supposing that the size of the 
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local averaging volume is equal to the cut-ofF scale Xk- 
A specific energy exchange rate is given by the ratio of 
the local energy exchange rate into a given volume to the 
local mass density. Using this definition we obtain the 
mean specific energy exchange rate which is as low as 
~ 10~^^nergs~^ cm~"^. However, its expectation value 
is proportional to the density. Thus, in regions as dense 
as the galactic halo the energy input due to the energy ex- 
change may be as large as « 10"^'' nergs"^ cm~^, which is 
comparable with the kinetic energy feedback by supernova 
averaged over the same galactic volume. 

On the basis of our results we consider a speculative 
picture: We assume that the injected energy propagates 
to MC scales and is therefore available for driving the MC 
turbulence leading to star formation. On the ground of 
this quite heuristic assumptions we propose an estimator 
for the star formation rate in a local volume in cosmolog- 
ical simulations: oc (ti'^/c^) to. 

The model reproduces approximately the observed global 
evolution of the star formation rate. Adopting this model 
means that merging and matter infall processes are con- 
sidered to enhance the star formation at least or to be the 
leading processes. 
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